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CALCULATION OF COEFFICIENTS OF THE OPTIMAL 
QUADRATURE FORMULAS IN THE W^''^~\0,1) SPACE 

Kh.M.Shadimetov, A.R.Hayotov 

Abstract. In this paper problem of construction of optimal quadrature formulas in VF2™'™~^^(0, 1) space is 
considered. Here by using Sobolev's algorithm when m = 1, 2 we find optimal coefficients of quadrature formulas 
of the form 

00 ; { 13=0 
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(3 . Many problems of science and engineering are reduced to a integral and a differential 

m 

equations or a system of such type equations. Solutions of such equations is often expressed 
<^ ■ with the help of defined integrals. But in many cases it is impossible to calculate these integrals 

exactly. Therefore approximate evaluation such integrals with possible high accuracy and less 
expenditure is one of the actual problem of computational mathematics. This is well demonstrated 
in example of theory of quadrature formulas. The task of finding of numerical value of one 
dimensional integral, in view of its well-known geometrical meaning, is often said the quadrature. 
There exist various methods of quadrature, which allow to calculate defined integral with the 
. help of finite number values of integrand. 

in ■ 

Present work also is devoted to one of such methods, i.e. to construction of optimal quadrature 
formulas for approximate calculation of defined integrals in the w!f^'"^ ^■*(0, 1) space. 
We consider following quadrature formula 



1 N 



(3=0 

with error functional 

N 

^Ar(x) = p(x)£[o,l](x) - '^CpSix - Xp) (2) 

13=0 

on space. Here Cp are the coefficients and Xfs are the nodes of the quadrature 

formula (1), £[0,1] (x) is the indicator of the interval [0, 1], p{x) is a weight function, S{x) is Dirac 
delta-function, if{x) is a element of the Hilbert space ^K— 1)^0,1), norm of a function in 
which is given by formula 

1/2 

Mx)\wt"'-'\0,1)\\ = { I {^^"^\x) + ^^^-'\x))' dx' 



Difference 

1 N 

{iN{x),(^{x)) = / p{x)ip{x)dx -^Cpip{xp) = £N{x)ip{x)dx. 

13=0 

is called by error of the quadrature formula (1). 

Quality of the quadrature formula (1) is estimated by norm of the error functional 

Wi^ixWt'^-'^iO,!)]] = sup \(e^(x),^(x))\. 

Mx)\w^"^'"^-'Ho,i)\\=i 

The norm of the error functional iN{x) depends on the coefficients Cp and the nodes x^. Choice 
of coefficients in fixed nodes is linear problem. Therefore we will fix the nodes Xj^ and the norm 
of the error functional will minimize by the coefficients, i.e. we will find 

II L '"^-^^*(0,l)|| =inf ||£^(a;)|iyf '"^-^)*(0,l)||. (3) 

C/3 

o o 

If it is found || In then said, that the functional £jv (x) correspond to the optimal 

quadrature formula in W2^'^~^\o, 1). Thus, we will get following problems. 

Problem 1. Find norm of the error functional £n{x) of the quadrature formula of the form 
(1) inWt'"'''^\QA) space. 

Problem 2. Find such values of the coefficients C^, which satisfy the equality (3). 

Problem 1 was solved in [1] and for square of the norm of the error functional (2) was 
obtained 

iKiv(x)ir = (-ir 



TV N 

|2 ' 



l3=0 i3'=o 

1 1 



/?=0 



(4) 



where 



e — e ^ \ 

Furthermore, the error functional (2), as shown in [1], satisfies following conditions 



(£jv(a;),x") = 0, a = 0,m-2, (6) 

(£A,(x),e-^) =0. (7) 
The norm of the error functional is multidimensional function of the coefficients (/3 = 



0, N). For finding point of condition minimum of the norm square of the error functional (2) in 



2 



conditions (6) and (7) in [1] applying the method of Lagrange undetermined factors following 
system was obtained 

; 

C^i^rnixp - X^) + Pm-2{Xi3) + de = / p{x)ll^rn{x - Xp)dx, f3^0,N, (8) 
n J 



7=0 



^CjX^ ^ p{x)x"dx, a^0,m-2, (9) 

n ^ 



7=0 





1 



C^e-^^ = / p{x)e-''dx, (10) 

7=0 { 

where ipmix) is determined by formula (5), Pm-2{x) is unknown polynomial of degree m — 2,d 
is a constant. 

In [1] also was proved existence and uniqueness of solution of the system (8)-(10). 

In present paper we will solve the system (8) - (10), i.e. we will solve the problem 2. 

Here interestingly that if system (8) - (10) is solved and the optimal coefficients Cp are found, 
then obtained quadrature formula will be exact not only for polynomials of degree m — 2, it 
also will be exact for function e~^. 

Direct solution of the system (8)-(10) is very hard. Therefore we will use Sobolev's algorithm 
[2], which is applied in solution such systems. 

Further, we briefly present this algorithm for the system (8) - (10). For this we will determine 
Cp for all integer values of /?, assume that for /3 < and {3 > N coefficients Cp are equal to 
zero. Let, furthermore, xp — hp, h — j^, N — 1,2, .... 

The system (8)-(10) by using the definition of convolution of discrete functions we reduce 
to the following form 

^m{hf3) * Cp + Pm-2m + de-'^^ = fm{hf3), hp e [0, 1], (11) 

= when ^ [0,1], (12) 

AT 

J2Cpihpr^q^, a = 0,m-2, (13) 

13=0 

N 

Y,Cpe-^^ = qm-u (14) 

/3=0 

where Pm-2{h/3) is unknown polynomial of degree m — 2, d is a constant, 

1 

fm{h(3) = p{x)ilJm{x - hl3)dx, (15) 







1 1 

Qa = f p{x)x°'dx, (a = 0, m — 2), Qm-i = J p{x)e~^dx. 



Consider following problem. 

Problem A. Find function Cp and polynomial Pm-2{hl3), unknown coefficient d, which 



satisfy the system (11) - (14) in given fm{h(3) and Qa {a — 0,m — 1). 

Idea of solution of the system (11) - (14), as in [2], consists on replace unknown function Cp. 
Namely, instead of (Tg we will introduce into consideration functions 



and 



u{h(3) = v{h(5) + Pm-2{h(5) + de-''^. 



In such substitution we need only express Cjj with the help of u{h(3). For this we must construct 
such operator Dm{h(3), which satisfies following equality 



Dm{hp)*il;mm^S{hp), 



(16) 



where S{h(3) is equal to zero when hp ^ and is equal to 1 when hp — 0. In works [3,4] discrete 
operator Djn{hp) was constructed, its properties were investigated and following theorems were 
proved 

Theorem 1. Discrete analogue of the differential operator — ^JL-i , which satisfies 
equality (16) have following view: 



Dm{h(5) 



(2m-2) 
P2m-2 



k=l 

m— 1 

-2e'^ + E A,, \P\ = 1; 

k=l 
m—1 

2^^+E#, /^ = o, 



(17) 



k=l 



where 



C = 1 + (2m - 2)6" + e"'* + 



h . j2m-2) 
h , „2h , ^ F2m-3 



(2m-2) ' 
P2m-2 



Ak = 



2(1 - \,f--^[X,{e^^ + 1) - e\Xl + 1)]^^'^ 



^kP2m-2{^k) 

^2m-2(A) = pf™-2)A^ = (1 - e^''){l - A)2"^-2 - 2(A(e2'^ + 1) - e^{\^ + l))x 



(18) 
(19) 



2m-2 



s=0 



X 



/i3(l - A) 



2m-6 



h{l - \f^-' + E^{\) + ... + 



/t^'»-3i?2^_4(A) 

(2m -3)! 



(20) 
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P2m-2 \ vim-z' ^'"^ Coefficients of the polynomial -P2m-2(A), is a root of the polynomial 
-P2m-2(A), |Ajk| < 1, Ek{X) is Euler polynomial [5]. 

Theorem 2. The discrete analogue Dm{h(3) of the differential operator — J^2m-2 when 
m = 1, 2, 3 satisfies following equalities 

1) D^{hp) * = 0; 

2) Dm{h(3) * e-^P = 0; 

3) Dm{hp) * {hpi)^ — 0, as n < 2m — 3, i.e. convolution of Dm{hpi) with polynomial of degree 
< 2m — 3 when m — 2,3 is equal to zero; 

4) Dm{hf3) * i(jm{hp) = S{h/3), where ■0^(^/5) is determined by formula (5). 
Further, taking account of (16), we get 

C^^Drr,ihP)*uihP). (21) 

So, if we find function u{h(3), then the optimal coefficients are found from equahty (21). 

Now we will find explicit form of the function u{h(3). From equality (11) we get, that u{h(5) — 
fm{hp) when hp e [0, 1]. Since C^^O when hp ^ [0, 1], then 

Cp = Drn{hp) * u{h(3) = when h(5 ^ [0, 1]. (22) 

We will calculate convolution v{h(3) = ijjm{h(3) * Cp when h[3 ^ [0, 1]. 
Taking into account (12), (13), (14), we have 

<m-{ J^""'' ' ' " "^"^ ' (23) 

^J^q^_,-De-^f^-Q('-^-^\hp)-Q^.,m, /3>N, 



where 



^qm-i + De-''^ + Q^^^-^\h(3) + Qrn^2{hP), < 0; 
-^''^^-^- (W^-^-"(-l)-- 



Q(^-3)(/i/3) = i 



r m-\-l 1 -r 

2 J ^2k-l ,,„.2k-l-a( i \a 



k=l a=0 ^ ' 



^ (2k-\-a)\-a\ 



m-1 2fc-l ^x2fe_i_a,/ . 



(24) 



e e i^_:^:i, e^(m°. (25) 

fc=[I!i±l] a=m— 1 ^ ' 7=0 

1 ^ 

^ = 4 E Ct^'^- (26) 



We denote 



7=0 



Ql;:i2{hf3)^Pm-2ihP) + Qm-2m, a' ^ d + D , (27) 



Q';^Uh(3)^Pm-2(hP)-QmMhP), a+^d-D. (28) 
Then, taking into account (11), (16), (23), (27), (28), we obtain following problem 
Problem B. Find solution of the equation 

Dm{h(3) * u{h(5) = when h(5 ^ [0, 1] 

which has view: 



u{hp) = < 



0<P<N; 



(29) 



(+) 



Here polynomials (5^-2 (^^/^)) Qm-2{hP) constants a , a"*" are unknowns. If we will find 
them, then from (27), (28) we have 

Pm-2m = liQl^UhP) + Q^r^UhP)), d = ^[a- + a+], 

Qm-2m = l{Ql;;UhP) - Q^^^UhP)), D = l[a- - a+]. 
Unknowns Qml2{h(3), (5^-2 (^/^)) found from conditions (22). Thus, problem B 

and respectively problem A are simultaneously solved. 
We introduce following designations 



A 



-h-y 



hj 



(30) 



a 



2e^ 
V 



2e'' 

P 



qm-1 - + h) + ^2(1 + h) + a+e-'^-i - /„(1 + h) 



(31) 
(32) 

(33) 



In present paper, using above mentioned algorithm [2], following theorems are proved. 
Theorem 3. The coefficients of the weight optimal quadrature formulas of the form (1) in 



the Wr 



(m,m— 1) 



(0, 1) space have view: 



^ m—l 

DmihP) * fm{hP)\p=0 + a + E + hXk), = 0; 

fe=l 

m—l 



(34) 



k=l 



m—l 



Dm{h/3)*fUhp)\/3=N + b+ EM + h), I3^N. 

k=l 



where Dm{h(3) is defined by (17), fm{hP) is defined by (15), ak, bk, a, b are determined by 
(30)-(33), Afe are roots of the polynomial -P2m-2(A), which given by equality (19) and \\k\ < 1. 

Theorem 4. The optimal coefficients of the quadrature formulas of the form (1) in the 
W^'^\{i, 1) space when p{x) — 1 have following view: 



Ca 



where h = N = 1, 2, .... 

Theorem 5. The optimal coefficients of the quadrature formulas of the form (1) in the 
W2 ' (0, 1) space when p{x) — 1 have following view: 



l-^-K(/i)(Ai-Af), /5 = 

h + K{h) ((e'^ - Ai)Af + (1 - Aie'^)Af-'^) , (3 = l,iV- 1 
-1 + e'^ - K{h){\, - Af )) , = TV. 



where 



K{h) 



{2e^ -h){\i-l) 



Ai 



2(e'^-l)2(Ai + Af+^) 

+ 1) - , 2'' + 1 - (( '' - l) v///-(( '' + l)-' + 2//(l-e'0 



1 _ e2/i + 2he^ 



, |Ai|<l, 



/i = ^,7V = l,2, 



Proofs of theorems 

Proof of theorem 3. We obtained, that 

Cp^Dm{h(3)*u{hl3). 

Hence, taking into account (17), (29) and properties of the operator Dm{hj3) we have 

00 

C/s = Dm{hl3) * u{hl3) = J2 Dm{hl3 - hj)u{hj) = 



7=— 00 



7=1 



-/17 



+ 



N 



+ J2 Dm{hp - h^)fm{h^) + Dm{h{^ + N) - h(3) 

7=0 7=1 

+QL'"22(l + /^7) + a+e-^-'" 



,^7+1 



"?m-l — Q 



(2m-3), 



1 + /i7) + 



7 



According to the definition of convolution of discrete functions we get 



—;—Qm-i+ 

7=— 00 7=1 

00 

+g(2--3)(_/,^) + Ql^l^(-h-f) + a-e^-^ - U-hl)] + J] D^{h{^ + N) - h(5) 



X 



7=1 



X 



ohl3 



Qm-i - Q('™-')(1 + hi) + QiJ^Uhl + 1) + a+e-'^^-i - U{hl + 1) 



7=1 



?i7 



7=1 

Hence when /3 = using equality (17) we obtain 



Co^D^{hP)*fUhP)\f3=o + 



-2e^ r e-^ 



P 



k=i ^''P U 



m— 1" 



4 

m— 1 

y 

^^^^ ^ 

_g(2m-3)(i + /i7) + Q(+)(1 + h^i) + a+e"^-'^^ - /^(l + /17)] . 
From (36), taking account of (30), (31) and (32) we have 



m—l 



Co = Drn{hl3) * fm{h(3)\p=o + a + ^(afc + b^)- 



k=l 



Let now 1 < /3 < iV — 1, then from (35) and (17), we get 



00 m— 1 



A 



7=1 ik=l P 



k \/3+7-l 



3-?i7 



00 m—l 



7=1 A;=l 



3/17 



-Qm-V 



(35) 



(36) 



Hence keeping in mind (30) and (31) we obtain 



m— 1 



k=l 

Now consider case when (3 — N . From (35), taking account of equality (17) and designations 
(30), (31), (33), we get 



oo m— 1 



A, 



7=1 k=l 



-h-y 



oo m— 1 



+Q^:U^ + h) + a^e-'-' - Ul + h)\+J2Yl Xt^^ 

7=1 fc=l 

_g(2»-3)(i + + g(+)^(i + /,^) + a+e-'-^-^ - U{1 + h^) 

m—l 



-Qm-1~ 



k=l 

Thus, we get the statement of theorem. Theorem 3 is proved. 

Consider case when p{x) = 1. Then system (11)-(14) of Wiener -Hopf type has following 
view 

r N 

7=0 

N 

EC7(/i7)" = ^, « = 0,m-2; (37) 

7=0 
N 

K 7=0 

where il^mix) is defined by equality (5), Pm-2{h(3) is polynomial of degree m — 2 of h(3, d is 
unknown coefficient 



Jmyiip) — ^ / ^ 

fc=l 



(38) 



Proof of theorem 4. Let m = 1, then from (37) we have 



N 



EC, 

7=0 



sign(/^/3- ^7) ,,fe/3 
4 



g/»/3-fe7 _ e/^7-/^/3) + ^e-/»/3 = /i(/i/3), /3 = 0, iV, 



^ (776"'^^ = 1 - e-^ 
7=0 



(39) 



(40) 



According to (34) in this case the optimal coefficients have following form 



C,= { D,{h(3)*h{h(3), /5 = l,iV-l; (41) 

D^{h(3)*h{hl3)\p=M + h, (5^N. 

To find optimal coefficients (41) we need to calculate the convolution D\{h[3) * f\{h[3). For 
this using equalities (17), (38) and theorem 2 we obtain 



D^{hf3)*h{hl3) = Dr{hl3)* 



1 = 



oo 

= -D^{hf3) * 1 = - J] A(/i7) = -2Di(/i) - A(0) = 

7=— oo 

^ _2 . i(-2e'^) - i • 2C = --(-4e'^ + 2(1 + e^^^)) = '^K^- 

p e" + 1 



Then from (41) we have 



P 





-1) 




hi 


2(e'^ 


-1) 


e'»- 


hi 


2(e'^ 


-1) 



13=1, N -I 



(42) 



Hence one can see, that for finding of coefficients we need to determine a and h. For this directly 



using (42), (40), from (39) taking account of /3 = 1, A?" — 1, we have 



7=0 



7=0 



7=0 



It 



C'y 6 ^^^^^ €■ 



.hfS-h-y 



hy-hl3 



v7=0 



7=0 



1 
2 



2(e'^ - 1) 



,/i/3-/i7 _ f,h-f-hl3\ _ 



7=0 



^ ^ e'^ + 1 V 1-e-'^ 1-e'* 
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1 



^ I 2(e^-l) 



+ e 



-hp 



1 - e 



h3+li 



a + be + 



2(e'^+^ - 1) 



a 2e'* 
2 ^ 6^^ + 1 



1 - e-i 



a 

"2 ^ e'^ + l 



1 ^ o + fee , e'^+i - 1 



2(e'* + l) 



1. 



Substituting obtained expression for S into (39) we get 



a e 
2 ^ e'* + 1 



4e 



+ e" 



2(e'» + l) 



-1 + de 



-hp 



,1 + e 



-1 



1-e 



1. 



Hence and from (40) we obtain 



a _|_ e" e— 1 _ e+1 

2 e''+l 4e ~ 4e 



2(e''+l) 

ae + 6 + %r 



l+e 
4 



e- 1 



Solving this system, we will find unknowns a, b, d, i.e. 



(43) 



eh + i' e^ + r 

Thus, substituting (43) into (42), we obtain the assertion of theorem. Theorem 4 is proved. 
Proof of theorem 5. In m = 2 from (34) we obtain 

D2{hp)*f2ihp)\p=o + a + ai + b,X^, /3 = 
D2{hp) * f2{hp) + aiAf + b^X^''^, 1 < (3 < N - 1 

D^m * f2my=N + b + a.Af + 61, P^N 

First we will calculate the convolution D2{h(3) * f2{h(3). Using equalities (17), (38) and 
theorem 2 we get 

.^hp^^-hp^^i-hp^^hp-i_^ 



(44) 



D2{h(5)*f2{h(5) = D2{h(5)* 

{h^f + {l-hpf 



-) = -D2m 



» = -lz..(.,)*(W. 



(45) 



From equality (20) we obtain, that polynomial 

P2(A) = A'(l - e^^ + 2he^) + {2{e^^ - 1) - 2/1(6''^ + 1))A + (1 - e^^ + 2he^) 



11 



has two roots Ai and A2, A1A2 = 1 and |Ai| < 1 



Ai + A2 



2h{e'"' + 1) - 2(e2'^ - 1) 



, p = 1 - e^^ + 2he 



Further, 



(1 _ e2/^ + 2he'') 

00 00 

7=— CX3 7=— 00 

00 ^ 

2D2{h)h^ + 2 E D2{h^){h^f = 2/i2(-(-2e'^ + 



7=2 



P 



°° A 9^2 



_ 2^ 

P 

2^ 



Ai 1 - Ai t^r^l - Ai 

2A? 



7=1 



i A jn2 



i=0 

Ai 



2^ 



Ail-Ai Vl-Ai (1-Ai)^ 



2^ 



_2^.^^iAi(l + A, 



Ai (l^Ai)3 



Ai 1- Ai(l-Ai)2 
_ 2^ 



(1 - A,)= 



Hence, taking into account equalities (19) and (46), we have 



2P 

p 



-2e' 



2(Ai)^(Ai(e^^ + l)-e^(A^ + l))(l + Ai) 
(A?-l)(l-Ai)3 



2^ 
P 



-2e 



, 2(Ai(e^^ + l)-e^(A^ + l)) 
(1-Ai)^ 



(1-Ai)^ 



4^ 
P 



- 2Aie^ + e^A^ + Ai(e^^ + 1) - e^(Af + 1) 

a^"Ao^ 



4/^2 (gh _ 1)2;^^ -4/12(6^^ - 1)2 



P (1-Ai)2 

4/12(6^^ - 1)2 1 

p A2 - 2 + Ai 

4/i2(e^ - 1)2 



p A2(l - 2Ai + Af) 
4/i2(e'^ - 1)2 1 



r, 2ft(e^fe+l)-2(e2h-l) _ „ 

^ l_e2h+2/ie'> ^ 



p 2/i(e2^ + 1) - 2(e2'* - 1) - 2(1 - 2e2'*) - Ahe^ 
Mi^{e^ - 1)2 4/i2(e'' - 1)2 



2/i(e2'^ - 2e^ + 1) 2h[e^ - Vf 



= -2/^. 
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So 

L>2(/i/3) * (hp)^ = -2h. (47) 
If we substitute equality (47) into (45) we obtain 

* hm = —D^m * {hpf = -^{-2h) = h. (48) 

Taking into account (48) from (44) we get 

h + a + ai + biXf, P = 0, 
Cp^ { h + aiXf + biX^-^, 1<P<N-1, (49) 
h + b + aiX^ + bi, p^N. 

Here ai,bi,a,b are unknown constants. If we find these constants, then we obtain explicit form 
of optimal coefficients. 

When m — 2 the system (8)-(10) of Wiener-Hopf type for optimal coefficients have following 
view 

^ sign(/,^-M (e'-^-'--e->->^ _ _ \ ^ ^ 

7=0 ^ ^ 

ehP + e-'*^ + e^-^^ + e^^'^ - 4 {hpf + (1 - /i/?)^ 



(50) 



4 4 

5^C, = 1, (51) 

7=0 

AT 

^C^e-^'' ^ I - e-\ (52) 

7=0 

To find explicit form of optimal coefficients, substituting equality (49) into (50), we obtain 
identity by variable hp. Equating to zero coefficients of e^^ , e~^^ , hp, (hp)^, we get system for 
unknowns ai,bi and solving it we will find ai,bi. 

First from orthogonality conditions (51), (52) coefficients Cq and Cjv we will express by C^, 



7 = 1,N -1 

N-l 

Co + Cjv = 1 — C-y, 



Hence 



7=1 

N-l 

Co + Cnc-^ = 1 - - ^ C^e-'^T. 

7=1 



e - 2 _ 1 - ei-'^^ 



Co = ^ + J:C,^-^, (53) 



7=1 



13 



c, 



1 „ e^-'^^ - e 



N 



Now consider first sum of the equality (50) and taking into account definition of sign(a;) we get 

7=0 ^ 

JL / phfS-h-y _ h-y-hl3 \ 

= E^7 ^ (^/3-/^7) - 



7=0 

TV 



7=0 ^ 
/phP_p-hp \ ^ / hp-h-Y _ h-Y-hp 

Co ( / - +E^7 — (/^/?- M ) - 

\ 7=0 7=0 7=0 7=0 

Co ( / - /^/? j + E — ^-^ - ) - 

1 / hp -hp ^ ^ 

"2 " e-^) - ^ E ^-re"^ -hP + Y. C,hj ] . (55) 
\ 7=0 7=0 



Using (49) - form of optimal coefficients C^, 7 = 1, — 1, we calculate following sums 

= ^ ^^^^ = ^ I' - + "^T^aT + ^^T^aT 

e / I - e^-i Ai - Af e'^-i , Af - Aie'^-i 
'/i^^ - + ai ' \ + ^'1 ' 



e - 1 V e'' - 1 - Ai \ie^ - 1 



Af-l 



•52 C^ — = r U— 7 — + ai , \ ^ 



e-1 e-lV e'^-l e'^-Ai 
7=1 ^ 

Af-Aie^-^ \ e / Ai - Af , , Af - Ai 

+ O1 r 7 ; 7 1 - /i + ai— h hi 



Aie'^-l y e-lV' ' ' ' 1-Ai ' Ai - 1 

13 



P /php-h-y _ ph-y-hP \ 
^3 = E^7( ^ {hP-h^)j^ 



14 



7=1 7=1 
7=1 7=1 

Hence using properties of geometric progression we get 

ehi^ fl- e-'^/^ Ai - Af+^e-'^^ , Af - Af-^e"'^''^ 

/ e'' - e'^+'^z^ Aie'^ - Af+'e'^^+'» Af e'^ - Xf-^e^+^f^ 
h— j— + ai^- + bi-^ ' 



\-e^ ' 1 - Xief^ ^ Ai - 

-h(3 ( /i/? + ai^— ^ + h \^_\ 1 + ^ n 



/5(Af+2 - A^^) + Ai - Af+^ , Af+i - Af+^-^ + /3(Af-'^ - Af+^-^) 

+ai/i t:; h Oi/l- 



2 



(l"Ai)2 (Ai-1) 
e^P [ h ^ aiXi ^ 6iAf A e^'^'^ / he^ ^ 



eh _ 1 -Xi Aie^ - 1 y 2 V 1 - 
Xie^ , Afe'* \ 1/ /i aiAf+^ 



l-Aie'^ Ai-e'^/ 2Vl-e^ Ai - e'* 
6iAf-'^ X 1 / /ie'^ aiAf+^e'^ feiAf-'^e^^ 



l-Aie'*/ 2ll-e'» 1 - Aie'' Ai - e'^ 
. « ^ Ai - Af +^ ^ Af - Af ^2 (/3 + 1)/? 

(T^AO^ + (vrTj2 

W-l N-l 
7=1 7=1 

e^-e Xie^ - Af e , Af - Aic 



S, = hY,C,^ = hY,{h + a,Xj + bX-'h = ^\ ' + 

7=1 7=1 

^ Ai - Af +^ - iVAf (1 - Ai) , ^ ^ Af +^ - Ai - jV(A^ - AQ 

P^T^T? + ' 
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Now by using equalities (53), (54), (55) and sums ^i, 5*2, S3, S4, S5, equating coefficients in front 
of e^^ and h(3 in left and right sides of equality (50), respectively, after some calculations we 
obtain following linear system for ai and bi. 

(e'^-l)(Af+^-Aie) , (e'--l)(Ai-Af +^e) _ 2-h h 
"l(l-Ai)(e-l)(e''-Ai) (e-l)(Ai-l)(Aie''-l) ~ 2 i?^' 

«i^ + &1aT^ = 0- 
Hence, solving system, we have 

_ (2e^-2-/.e^-/.)(A,-l) , 

- 2(e'^-l)2(A, + Ar^) " ^^^^ 

_ (2e^^2-/^,e^-M(Ai-l) , 

2(e'^-l)2(A, + Ar^) ^1-^^^)- 

From (53), taking account of Si and ai,bi for Co we get 



e 



(2e''-2-/ie'*-/i)(A2 + Af -Ai-Af+^) 



2(6^^ - l)2(Ai + Af +^) ■ ^^^^ 



And from equality (54), taking into account 5*2, (56), (57) we obtain optimal coefficient C, 



N 



_ he^-e^ + l e^(2e^ - 2 - he^' - h) (A^ + Af - Ai - Af +^) 

e'^-l 2(e'^-l)^(Ai + Af+i) " ^^^^ 

We denote by 

_ {2e^-2-he^-h){\,-l) 
^ ' 2(e^ - l)2(Ai + Af +1) 
and taking into account equalities (56), (57), (58) and (59) we obtain assertion of theorem. 
Theorem 5 is proved. 
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